Spin Excitations in La 2 Cu04 : Consistent Description by Inclusion of Ring-Exchange 
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We consider the square lattice Heisenberg antiferromagnet with plaquette ring exchange and a 
finite interlayer coupling leading to a consistent description of the spin-wave excitation spectrum in 
La2CuC>4. The values of the in- plane exchange parameters, including ring-exchange Jn , are obtained 
consistently by an accurate fit to the experimentally observed in-plane spin-wave dispersion, while 
the out-of-plane exchange interaction is found from the temperature dependence of the sublattice 
magnetization at low temperatures. The fitted exchange interactions J = 151.9meV and Jn = 
0.24J give values for the spin stiffness and the Neel temperature in excellent agreement with the 
experimental data. 
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The magnetic properties of La2CuC>4 have been the 
subject of many detailed investigations over the last 
decade. Understanding this undoped parent compound 
of high temperature superconducting cuprates is a pre- 
condition for the many theories which describe metal- 
lic cuprates by doping carriers into a layered antiferro- 
magnet. The conventional starting point for undoped 
cuprates is the two-dimensional (2D) spin-1/2 Heisen- 
berg model with only the nearest-neighbor exchange in- 
teraction J jy], which thereby provides the important 
magnetic energy scale needed as input to theories for 
the metallic and superconducting properties of doped 
cuprates. Despite the substantial progress on the the- 
ory of the 2D Heisenberg antiferromagnet B, which in- 
cludes such physical properties as the temperature de- 
pendence of the magnetic correlation length [|,||, some of 
the experimental facts for La2Cu04 have clearly demon- 
strated that a complete description of the magnetic ex- 
citations requires additional physics not contained in the 
2D Heisenberg model with J only. Examples include the 
asymmetric lineshape of the two-magnon Raman inten- 
sity |j| or the infrared optical absorption 0], which have 
led to proposals that spin-phonon interactions @, res- 
onant phenomena purely fermionic contributions 
fiof , or cyclic ring-exchange 11 -13|| need to be included. 
In particular, the importance of ring (plaquette) ex- 
change has very recently found direct experimental sup- 
port from the observed dispersion of the spin- waves along 
the magnetic Brillouin zone boundary |l4[ . 

The 4-spin plaquette ring-exchange interaction J n was 
considered rather early as a possible non-negligible cor- 
rection to the nearest-neighbor Heisenberg model Jl5],[l6| . 
This higher-order spin coupling arises naturally in a 
strong coupling t/U expansion to fourth order for the sin- 
gle band Hubbard model at half-filling Its quan- 
titative significance was recently demonstrated in high 
order perturbation expansions and ab initio cluster 
calculations in realistic three-band Hubbard models for 
the CuC>2 planes j2C]] . These derivations of effective spin 
models for the low-energy magnetic properties of the un- 



doped CuC>2 planes led to the estimate Jn/J = 0.11. 
A linear spin- wave analysis of the spectrum in La2Cu04 
at 10K in Ref. |14| has deduced a considerably larger 
value Juj J — 0.41. The necessity for a sizeable Jn has 
recently also been conjectured for the spin ladder com- 
pound La 2 Ca8Cu2404i |2l[ . 

The spin- wave theory and the quasiclassical phase di- 
agram of the frustrated Heisenberg model with ring ex- 
change was investigated in Ref. Q . However, the quan- 
tum and thermal renormalizations have not so far been 
taken into account in the spin-wave theory. It is well 
known for the quasi-2D Heisenberg model without the 
ring-exchange term (see Ref. [^3| and references therein) 
that such renormalizations can substantially change the 
excitation spectrum of the system. The authors of Ref. 
|p^| considered the simplest renormalization of the spec- 
trum by allowing for an overall quantum renormalization 
factor, which was obtained for the 2D Heisenberg model 
with nearest-neighbor exchange within a 1/5 expansion 
to order 1/S 2 |24| and by series expansion from the Ising 
limit p5|. However, in the absence of a consistent anal- 
ysis of the spin-wave renormalization in the presence of 
ring-exchange or a next-nearest neighbor coupling, it is 
not possible to provide an accurate determination of ex- 
change parameters. Indeed, as we show below, the effects 
of quantum and thermal fluctuations are not simply cap- 
tured by a single renormalization factor, and a consistent 
treatment of the spin-wave spectrum with J n 7^0 to order 
1/S reveals that the previous early estimate J=136meV 
from high energy neutron scattering [ p6| or two-magnon 
Raman scattering |27| requires a correction at least as 
large as 10%. Also, we show that the recent estimate 
Jo = 0.41J in Ref. |Q appears to be twice as large as 
the value calculated by accounting systematically for 1/S 
renormalizations. 

In this paper we consider the corrections to the spin- 
wave spectrum to first order in l/S for finite J n using 
a self-consistent spin- wave theory |^3| . We obtain values 
of the in-plane and interplane exchange interactions of 
La2CuC>4 allowing an accurate fit of the dispersion. We 
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verify that the obtained exchange interactions correctly 
reproduce the measured values for the spin stiffness and 
the Neel temperature. 

We start from the Heisenberg model with ring- 
exchange pjniJii 



j j' j" 
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where J, J', and J" are the first (8), second (8') and third 
(<5")-nearest-neighbor in-plane exchanges, <5j_ connects to 
the nearest-neighbor sites in the adjacent planes with in- 
terplane exchange Jj_, and (ijkl) denotes the four sites of 
a planar plaquette involved in the ring exchange. We use 
the Dyson-Maleev representation for the spin operators 
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where A and B denote the magnetic sublattices; a\,cii, 
and bj,bj are Bose operators. After substituting Eqs. 
(||) and (j^j into the Hamiltonian (|l|), we decouple quar- 
tic terms into quadratic ones according to the procedure 
described in Ref. ]23|j . Keeping consistently all terms to 
order 1/5, we obtain 



H = S 
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Diagonalization of this Hamiltonian yields the spin- wave 
spectrum 



E± = ^Al-Bl, (7) 
A k = 4S[J 7 - J' 7 '(l -4)- JuS 2 (^ + 7 K')] 



-AJ"Si'{\ -vi) + 2J ± S 1± _ , 



(8) 



B k = 4S(J 7 - J n S 2 7s D K + 2J x 57j.i^- l , (9) 
where 

= ( cos + cos ■ = cos ^ cos > 

= (cos 2*^ + cos2fc y )/2 , = cosfc 2 , (10) 

and 7 = 7 5, 7 ' = 75' etc.; the lattice constants are set 
to unity. Since the equality = 2 7 ^ — 7^ is satisfied, 
the spectrum given by Eq. u7|) is necessarily gapless. It 
is apparent from this result for the dispersion that the 
renormalization coefficients { 7( j, 7 ^} cannot be combined 
into a single overall renormalization factor. The averages 
of the bosonic operators which enter in Eqs. (||) and (|^) 
are 
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The expression for the sublattice magnetization reads 

S = S-{a\ ai ). (12) 

As previously discussed for quasi-2D magnets pjj , Eqs. 
@i @> an( l (0) mus t be solved self-consistently. The 
spin- wave velocity c is obtained by expanding the disper- 



sion at small wavevector k = y fc| + ky, leading to 

E k ~ck, c= 2V2Z C JS , (13) 

where 



where Di = a% for % G A and Di = b\ for i £ B; also we 
use the notation Js = J, J51 = J' etc. The renormaliza- 
tion of the bare exchange parameters due to quantum- 
and thermal fluctuations is described by the coefficients: 



1 - [(ajai) + (aib i+d )]/S, d = 8,S± 
1 - [(aj ai ) - (ata i+d )]/S, d = 8',8" 
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is the spin-wave velocity renormalization factor. For the 
spin stiffness we obtain the result 



p s = JS 2 Z p , 

Z p = (AJSo/A )Z 2 



(15) 
(16) 



Given c and p s , the transverse susceptibility follows im- 
mediately as 

X± = Ps/c 2 = So/(2SA ) = Z X /(8J) (17) 



2 



where Z x = Z p /Zl. For J' = J" = J n = J ± = 0, i.e. 
the 2D Heisenberg model with only nearest-neighbor ex- 
change, the self-consistent numerical solution of Eqs. (||) 
and (jil|) at zero temperature gives |^ , ^9| , ^3| 7 = 1.158 
and Sq = 0.303, which corresponds to Z c = 1.158, Z p = 
0.702, and Z x = 0.524. We note that due to the re- 
lations ([l6]) and (|l7j), Z p and Z x contain partially con- 
tributions of order 1/S 2 . The above numbers are close 
to those found for the 2D nearest-neighbor Heisenberg 
model in a systematic 1/S expansion to order 1/S 2 : 
Z c = 1.179, Z p = 0.724, and Z x = 0.514 |§. However, 
we show below, the presence of next-nearest neighbor and 
plaquette ring-exchange terms substantially alters these 
numbers. 

We use Eqs. (||), and ( [li"|) at zero temperature to 
fit the experimentally determined planar spin-wave dis- 
persion using the data at T=10K from Ref. H). To 
restrict the number of fitting parameters, we suppose 
J' = J". This restriction is well justified by the pertur- 
bation expansions of the half-filled one- and three-band 
Hubbard models Jl7| |l9t] - The inelastic neutron scatter- 
ing data together with our fit result along a selected path 
in the Brillouin zone at T = 10K are shown in Fig. 1. 
The best fit is obtained for the following parameter set: 

J = 151.9 mcV, J' = J" = 0.025 J, J n = 0.24 J. (18) 

While the values of J' and J" are practically indistin- 
guishable from those in Ref. Jl4j and J in (18) is only 3% 
larger, our extracted value of </□ is 50% lower. 

For the corresponding groundstate sublattice mag- 
netization Sq and for the renormalization parameters 
{idild I we obtain 

So = 0.319, 7 = 1.158, i = 0.909, 7" = 0.852 , (19) 
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2.220, 75 = 1.971, 7°, = 1.721 . 



In our notation, the spectrum used in Ref. jL4j corre- 
sponds to equal renormalization factors: 7 = 7' = 7" = 
7° ~ 1.18; the parameter values obtained with this spec- 
trum were J' /J = J" /J = 0.020, and J a /J = 0.41. 
However, some of the 7-coefficients show a remarkable de- 
viation from the value 1.18 found for the 2D system with 
nearest- neighbor exchange only Jm| . In particular the 
renormalization coefficients {7^} for the ring-exchange 
deviate very strongly. We emphasize again, although 
some fitting parameters are very similar to |14|, it is the 
self-consitently renormalized parameters {7^,7^} which 
allow us to obtain an accurate and reliable set of bare 
superexchange couplings. 

In the self-consistent spin-wave theory presented here 
the in-plane magnon spectrum varies only weakly with 
temperature at T <C J. Although the spectrum changes 
qualitatively in the same way as found experimentally, it 
accounts only for a few per cent of the observed changes 
in the zone boundary dispersion of the data at T = 295K 
in Ref. 1T1. 
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FIG. 1. Spin- wave dispersion along high symmetry direc- 
tions in the 2D Brillouin zone. The triangles are the exper- 
imental results of Ref. [14] for La2CuC>4 at 10K. The solid 
line is the result of a fit to the spin-wave dispersion result (7) 
leading to the exchange couplings as listed in Eq. (18). The 
dashed line is the fit of Ref. [14] . 



From the parameter values obtained above we deduce 
the spin stiffness, the spin-wave velocity, and the trans- 
verse susceptibility: p s — 23.8 meV, c — 206 meV and 
X± = 4.8 x lO^KT 1 . The corresponding values of the 
renormalization factors as calculated from Eqs. (pl|), 
dT|), and @ are Z c = 0.96, Z p = 0.63 and Z x = 0.68. 
These values differ substantially from those for the 2D 
nearest-neighbor Heisenberg model. We note that our 
value for the spin stiffness is in very good agreement with 
the earlier estimate p s = 23.9meV found from fitting 
the spin-spin correlation length £(T) at T > Tjv to the 
non-linear cr-model result Q with the correct preexpo- 
nential factor ||. This agreement strongly supports the 
validity of the self-consistent renormalized spin- wave the- 
ory. 

As discussed in Refs. [ p3pl| ], it is difficult to fit the 
value Ji from measurements of the out-of-plane spin- 
wave spectrum. Instead, we use an alternative proce- 
dure and fit the temperature dependence of the sublat- 
tice magnetization at temperatures T < Try/2 where 
the above theory is reliable (cf. Ref. |23| and refer- 
ences therein) employing the exchange parameters listed 
in (|l8|). In this way we obtain J±/J = 1.0 x 10~ 3 and 
7^/7 = 5.67 x 10~ 4 which practically coincides with the 
previous estimate in Refs. p3| , pi[ . 

As another test of the above results, we also calcu- 
late the Neel temperature for the obtained exchange pa- 
rameter values. On the basis of a renormalization-group 
approach and a 1/N expansion in the O(N) quantum 
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nonlinear cr-model, the result for the Neel temperature 
of a quasi-2D isotropic Heisenberg antiferromagnet has 
the form 



N 



47T/9 S 



T 2 

In 4^ 



3 In 



Ti 



0.0660 



(20) 



where a r = (7^/7)t=0j and c and p s are the respective 
groundstate spin-wave velocity and spin stiffness given 
by Eqs. (15) and (|T^). With the above parameter values 
( |l8| ) we obtain Tjv = 328K, in almost perfect agreement 
with the experimental value Tjv = 325X | Ej9 |. 

In conclusion, we have considered the renormalization 
of the spin- wave spectrum to order 1/S for the Heisen- 
berg antiferromagnet in the presence of plaquette ring- 
exchange. The results allow for an accurate fit of the 
magnon dispersion in La2CuC>4 and a consistent deter- 
mination of the exchange coupling parameters for this 
material. As an independent check of the parameter 
set, the spin stiffness and the Neel temperature are cor- 
rectly reproduced. With J = 151.9meV the obtained 
value for the bare ring-exchange coupling J n /J = 0.24 
in La2CuC>4 is significant although due to its 4-spin cou- 
pling the ring-exchange term in the Hamiltonian will be 
reduced by a factor S 2 with respect to the J term. The 
magnitude of J n suggests that for hole doped cuprates 
ring-exchange might be relevant, too, and we may postu- 
late that it is connected to recent proposals of staggered 
circulating currents in underdoped materials [ j32"|]33| . The 
role of ring-exchange for the spectral lineshape of the B\ g 
shift in Raman experiments and for infrared absorption 
remains to be reexplored and contrasted to the recent 
proposal of a purely fermionic origin of spectral weight 
at higher energies |Hj. 
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